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Abstract 

We compute supersymmetric QCD and supersymmetric electroweak corrections to the 
b and t quark masses and to the r lepton mass in the presence of CP phases valid for 
all tan (3. The analysis includes one loop diagrams arising from the exchange of gluinos, 
charginos and neutralinos. We find that the CP effects can change both the magnitude as 
well as the sign of the supersymmetric loop correction. In the region of parameter space 
studied we found that the numerical size of the correction to the b quark mass can be as 
large as 30 percent and to the r lepton mass as much as 5 percent. For the top quark 
mass the correction is typically less than a percent. These corrections are of importance 
in unified models since b — r and b — t — r unification is strongly affected by the size 
of these corrections. More generally these results will be of importance in analyses of 
quark-lepton textures with the inclusion of CP phases. Further, the analysis presented 
here is relevant in a variety of low energy phenomena where supersymmetric QCD and 
supersymmetric electroweak corrections on b mass enter prominently in B physics, e.g., 
in processes such as B® d — > l + l~ . 
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1 Introduction 



The b quark mass plays a very important role in theories of fundamental interactions as 
it enters in tests of b — r unification [1] and in b — t — r unification [2]. [For recent works 
on Yukawa unification in SU(5) and in SO(10), see refs.[3]]. In addition the b quark mass 
enters in a variety of low energy phenomena. The physical mass, or the pole mass M b , of 
the b quark M b is related to the running mass m b (M b ) by inclusion of QCD corrections 
and at the two loop level one has[l] 
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M b = (1 + + 12A^^L)m b (M b ) (1) 

Now m b (M b ) is derived from m b [M z ) by the running of the renormalization group equa- 
tions and we focus here on the quantity m b (M z ) which is the running b quark mass at 
the Z scale. This quantity can be written in the form 

m b (M z ) = h(M z )^= cos/3(l + A b ) (2) 

Here h b (M z ) is the Yukawa coupling for the b quark at the scale M Z) (3 is defined so that 
tan/3 =< H 2 > / < Hi > where H 2 is the Higgs field that gives mass to the up quark 
and Hi is the Higgs field that gives mass to the down quark and the lepton, and A b is 
loop correction to m b . In previous analyses supersymmetric QCD and supersymmetric 
electroweak loop corrections to the running b quark mass at the Z boson mass scale have 
been computed and it is found that in the large tan/3 region these corrections can be 
rather large[4, 5, 6, 7, 8]. In this paper we investigate the effects of CP phases SUSY 
QCD and SUSY electroweak loop corrections on the b quark mass. We then extend the 
analysis to include the effects of CP phases on the loop corrections to the top quark mass 
and on the loop corrections on the r lepton mass. It is well known that supersymmetric 
theories contain CP phases via the soft SUSY breaking parameters which are in general 
complex and thus introduce new sources of CP violation in the theory above and beyond 
what is present in the standard model. Thus, for example, in mSUGRA[9] the soft SUSY 
breaking is characterized by the parameters m , mi, A and B , where m is the universal 
scalar mass, mi is the universal gaugino mass, A is the universal trilinear coupling and 
Bo is the bilinear coupling. In addition mSUGRA contains a parameter // which is the 
co-efficient of the Higgs mixing term in the superpotential. After spontaneous breaking of 
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the electroweak symmetry it is convenient to replace Bq by tan f3 and the thus mSUGRA 
at low energy can be described by the parameters mo, mi, Aq and tan/?. In the presence 
of CP violation mSUGRA contains two CP violating phases which cannot be removed by 
field redefinitions. These can be chosen to be the phases of fj, and A , i.e., // = |/xo|e 26 * M 
and Aq = \Ao\e taA . In the more general scenario of the minimal supersymmetric standard 
model (MSSM) one has many more phases. In this analysis we will consider this more 
general situation of MSSM where in general we will allow for independent phases for the 
gaugino masses rhi so that 

rhi = \rhi\e^% % = 1,2,3 (3) 

where i=3,2,l refer to the SU(3)c, SU(2) L and U(l) gauge sectors. One problem that 
must be taken into account with inclusion of CP phases is that of the current experimental 
constraints on the electric dipole moments (edms) of the electron, of the neutron and of 
the Hg" atom. Thus the current limit on the electron edm is d e < 4.3 x 10~ 27 ecm[10] 
and on the neutron is d n < 6.5 x 10~ 26 ecm[ll]). Similarly the limit on the Hg" atom is 
also very stringent, i.e., dn < 9 x 10~ 28 ecm[12]. Typically, with phases 0(1) the SUSY 
contributions to the electron and the neutron edms are already in excess of the current 
experimental limits. There are many suggestions on how to overcome this problem. These 
include fine tuning[13], large sparticle masses[14], suppression in the context of Left-Right 
models[15], the cancellation mechanism[16] and suppression by the choice of phases only 
in the third generation [17]. Specifically in Refs.[16] it is shown that in a large class of 
SUSY, string and brane models one can have cancellations and the satisfaction of the 
edm constraints for the electron and the neutron. In Ref.[18] it was also shown that if 
cancellations occur at a point in the parameter space then this point can be promoted 
into a trajectory where the cancellations occur by using a scaling. Further, in Refs.[19, 20] 
the implementation of the cancellation mechanism including the H X99 atomic constraint 
has also been carried out. In some of the above scenarios it is then possible to have large 
phases and also satisfy the edm constraints. It is now well known that if the phases are 
large they affect a variety of low energy phenomena. Some recent works in this direction 
have included the effects of CP phases on the neutral Higgs boson system [21, 22, 23], on 
collider physics[24, 25, 26], on g-2[27], in B° sd -> l+l~[20, 28] and in other B physics[29]. 
The number of such phenomena is rather large and a more complete list can be found in 
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Ref. [30]. In the analysis below we will assume that the phases are indeed large and that 
the satisfaction of the edm constraints occurs because of one of the mechanisms discussed 
above. 

The outline of the rest of the paper is as follows: In Sec. 2 we discuss the procedure 
for the computation of the SUSY QCD and SUSY electroweak correction to the fermion 
masses and compute explicitly the corrections to the b quark mass arising from the gluino, 
chargino and neutralino exchange loops. The analysis given is valid not just for large tan (3 
but for all tan/3. In Sec. 2 we also discuss the limit when tan/3 is large and the phases 
vanish and show that our analysis in this limit is in agreement with the previously derived 
results. In Sec. 3 we discuss the SUSY QCD and SUSY electroweak corrections to the top 
quark mass with the inclusion of phases arising from the gluino, chargino and neutralino 
exchange loops. In Sec. 4 we give the analysis for the SUSY electroweak correction to the 
r lepton mass arising from the exchange of charginos and neutralinos. As in the case of 
the analysis of the correction to the b quark mass, the corrections to the top mass and the 
r lepton mass include all allowed CP phases and the analysis is valid not just for all tan (3. 
The analysis given in Sees. 2, 3, and 4 is the first complete one loop analysis of the SUSY 
QCD and SUSY electroweak corrections to the b,t and r lepton masses with inclusion of 
phases and without any apriori limitation on tan/3. The analytic analysis of Sees. 2, 3 and 
4 contain the main new results of this paper. In Sec. 5 we give a numerical analysis of the 
size of the loop corrections to the b, t and r lepton masses. Sec. 6 is devoted to summary 
and conclusions. 



2 Correction to the b Quark Mass 

We begin by exhibiting the general technique we use in the computation of the corrections 
to the b quark mass. Our procedure here is similar to that of Ref [7, 8]. In this procedure 
one first evaluates the loop corrections to the Yukawa couplings to the Higgs boson. At 
the tree level the b quark couples to the neutral component of the Hi Higgs boson (i.e., 
Hf) while the couplings to the H 2 Higgs boson is absent. Loop corrections produce a shift 
in the H® coupling as expected and in addition also generate a non-vanishing effective 
coupling with H®. Thus the b quark coupling with the Higgs in the presence of loop 
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Figure 1: Exhibition of a generic supersymmetric loop contribution to the b quark mass. 
All particles in the loop are heavy supersymmetric partners with S±, S 2 being heavy scalars 
and / a heavy fermion. 

corrections can be written as [7, 8] 

-L bbH0 = (h b + 5h b )b R b L H° 1 + Ah b b R b L H° 2 * + H.c. (4) 

where we are using the normalization on the Higgs fields of Ref.[31]. The correction to 
the b quark mass is then given by 

Re(Ah b ) Re5h b 

A b =[ tan/5+— — ] (5) 

n b n b 

We now return to the details of the analysis. In Fig. 1 we exhibit the typical supersym- 
metric loop correction, which is either a supersymmetric QCD correction or a supersym- 
metric electroweak correction, that contributes to the b quark mass. The basic integral 
that enters in the computation of Ah b and 5h b is 



f d 4 k 7TIJ+ ft 

J (2tt) 4 {-m) + k 2 ){-m\ + k 2 ){-m\ + P) 



In the approximation that the external momentum can be set to zero the ]k term in the 
numerator inside the integral in Eq.(6) can be neglected and in this case the integral is 
given by 



(4vr) 2 f Sl 52 
where the function / (m 2 , mf , m 2 ) is given by 
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m~ 

j{m\m\.m\\ (7) 
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ft 2 2 2\ ^ / 2 2; m j , 2 2j m , 2 2r m i \ 

J (m 2 — m 2 )(m 2 — m 2 )(m 2 — m 2 ) J m 2 mf J m 2 - 

for the case i 7^ j and 

/(m 2 , m 2 , m 2 ) = j ^ 2)2 (mHrJ^ + (m 2 - m 2 )) (9) 

for the case i=j. In general one will have gluino, chargino and neutralino exchanges in 
the loops. We compute the contribution from each of these below. 

2.1 Gluino Exchange Contribution 

Fig. 2 exhibits the relevant loop involving g and &i )2 exchanges. The gluino interaction is 
given by 



-L S = V2g s ±± T^(-b"^g a e-^% + b^-^~g a e^W L ) + H.c. (10) 

where a, (3 are the quark color indices taking on values 1,2,3 and 

2 _ _ 2 
b L = D biibi, b R = J2 D b2i bi (11) 



i=i i=i 



Here D^j is the matrix that diagonalizes the b squark mass matrix and bi are the b squark 
mass eigen states. For the computation of Ahj, we need the bbH® interaction which is 
given by 
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^ = EE G «M+ ff - & ( 12 ) 

i=i i=i 

Using the interactions of Eqs.(10) and (12) and the identity 



one finds 



a a 1 1 
T ap T "s = 7, Pe<7 <W - o <W M ( 1 3 ) 

a=l Z 



A^ = -EE ^V^m,G^^/(m 2 , m? , m? ) (14) 
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Figure 2: One loop contribution to the bottom quark mass involving exchange of gluino, 
charginos and neutralinos in the loop. 



where 



<jXfz {(-l + l^n 2 9 w )D; ii D blj -lsm 2 9 w D* b2i D b2j }smP 



V2 2cos^ LV 2 ' 3 3 



+ 2M^cos/A^ (15) 



For the computation of Sh b one needs the bbH® interaction which is given by 

2 2 

W° = E E HifbfbjH* + H.c. (16) 

i=l j=l 



From Eqs.(lO) and (16) one finds 



2 2 2 

^ = - E E -^e-^m^D^D^fiml m? , m? ) (17) 
8 =j j=i 



where 



= + i Sin2 M^m^ " ^ sin2 °wD* m D b2j } cos/5 

-2M w cos(3 [DmDbl > + " 2M w cosf3 DlaiDbli (18) 



2.2 Chargino Exchange Contribution 

Fig. 2 exhibits the relevant loop involving x~ and exchanges. By carrying a similar 
analysis of the gluino exchange, the chargino contribution gives for Ah b the result 

2 2 2 m , 

A ^ = " E E E g'EUVCiDh - hV;,D^{k b Ul,D nj )^ 2 f{ml p m|, m\) 
i=i j=i fc=i 
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2 2 2 m v +m v + 
-EEEWW - W 2 ^ 2 J(^^ 2 Ai fc )^-^/K,m^,m^ + ) (19) 

j=l jr' = l fc=l 



where 



^2 = " I Sin2 W D ™ D * + I Sin2 ^^A 2j } sin/5 



—- [D* tli D tlj + D* t2i D t2j ] - » J ' l - D* t2i D tlj (20) 



and 



and where 



V2- 2sk?^- Q ^- R ^ < 21 > 



Qij = y 2^2^-1 

Rij = ^[^^i^-i + /i*f/,2^ 2 ] (22) 



and 



*«> = n '""'I m (23) 

y Imw sin p (cos p) 

The matrices U and V diagonalize the chargino mass matrix by a biunitary transforma- 
tion. Similarly the matrix D ti j diagonalizes the stop mass matrix. 



The computation of 5hb gives 

2 2 2 m 

16tt 2J V ""^ 



i=l j=l k=l 

E E E - w^}(fc6^ 2 Ai fc )^^^/K. < + ) ( 24 ) 

1=1 j=l fc=l 



where 



F i7 - oM z r/ l 2 , N „ 2 , 



{(o ~~ q sin °w)D* tli D tlj + - sin 9 w D* t2i D t2j } cos/3 



V2 2cos#vk 2 3 ^-u*^ 1 3 



+ 2M w smp DtllDt2 > (25) 
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and 



T2 = -l Q " (26) 

2.3 Neutralino Exchange Contribution 

Fig. 2 exhibits the relevant loop involving neutalinos (x k , A: = 1 — 4) and sbottom (61,2) 
exchanges. The contribution of these loops to Ah b is 

2 2 4 

' ■>('.'. \(\, < n. - ^..d^axr* nt. -1- ^..n^.l 

IGtt 2 -' v x *' 



M ^-EEE 2G?.{a 6fc £> Hj - - l bk D b2j }{(3* bk D* bli + a bk D* m }^f(m 2 x o, m? , m|) 
i=i j=i fc=i 



4 4 2 m r om y o 

- E E E 2r y {ay£> 61Jfc - TfcAad^u* + ««^w*}-S^/K' m x°' ( 27 ) 
i=ij=ifc=i i07r 1 J 



where 



gm b X 3k 

a b k 



2mw cos (5 



(3 bk = eQ b X'* k + -JL^x'* k (T m - Q b sin 2 

COS U\\r 

-f bk = eQ b X lk — — X 2k (28) 

COS t>w 



and where 



X' lk = X lk cos 9 W + X 2k sin 9 W 

X' 2k = -X lk sin 6 W + X 2 k cos 6 W (29) 

Here T 3b = —1/2, = —1/3 and e = gsin#ty and the matrix X diagonalizes the 

neutralino mass matrix (see, e.g., Ref. [32] for notation and definitions). appearing in 
Eq.(27) is defined by 

^l=-2ib [ ^- Q « cos " (30) 

where 

^ = 2lb KXl * A * j + *t x « x v - + ( 31 ) 
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Finally the analysis of the same loops gives the neutralino exchange contribution to Shb 

2 2 4 mo 
* = "EEE 2H ij {a hk D hlj - lhk D h2j }{l3l k D* m + a hk Dl 2 ^-^f{m\^m\,m\) 

i=l j=l fe=l iD7r 

4 4 2 m om o 

-EEE 2Ay{a w £> 6 i* " TkAakH/^i* + »MD* b2 J^-^/K fe , mjo, m^ ) (32) 
i=ij=ifc=i iD7r J 

where 

A similar analysis for the corrections to the top quark mass is given in Sec. 3 and for the 
r lepton mass is given in Sec. 4. 



2.4 The Limit of Vanishing Phases and Large tan/3 

We check the results of Sec. 2. 1-2 with previous analyses by taking the limit when phases 
vanish and when tan (3 is large. We first look at the gluino exchange contribution. Here 
in the large tan (3 limit in Eq. (5) the first term dominates. Further, setting £3 = and 
using the approximation of ignoring the sbottom mixing we find that the correction to 
the b quark mass from the gluino exchange from Eq.(14) is given by 

Af = ^itan/3/K,m| 2 ,M|) (34) 

This is precisely the result obtained in previous analyses for the case with large tan (3 and 
no phase (see Eq. 6 of [4]). Next we consider the chargino exchange contribution given by 
Eq.(19). To compare the results with previous analyses we assume that \fi\ » Mw,fh2- 
In this limit the matrices U and V are unity, and only the terms i=2, j=l, k=2 in the 
sum contributes. Further, using the approximation that the matrix D is unity, putting 
m + = u and setting the phases to zero one finds 

Am f = £ t&n p E2lkbk » f ^ :Tn 2 m |) (35) 
rib Lbit 3 

From the definitions of kb, k t , and in the limit of large tan (3 and having the relation that 

m A t » M z one finds that our chargino result with no phases gives 

ax+ _ 9 2 m 2 t m A t y 2 2 2 

A& " 327r 2 M^sin/3cos/5 / ^' mt2 ' /i } W 
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This is the result derived in previous analyses in the large tan/3 limit in the absence of 
phases(see Eq. 7 of [4]). In both of the above cases, i.e., for the gluino exchange and for 
the chargino exchange we find that the loop correction to the b quark mass increases with 
tan (3 and thus it can get large for large tan (5. 



3 Correction to the top Quark Mass 

In this section we discuss the loop corrections to the top quark mass in the presence of 
CP phases. The analysis proceeds in a manner very similar to the analysis of the bottom 
mass correction. Thus the running of the top mass at the Z scale is given by 

m t (M z ) = h t (M z )^=sm(3(l + A t ) (37) 
where A t gives the loop correction to m t . The loop correction A t is then given by 

At=i ^ihi catp+ ^m ) ,38) 

We carried out a one loop computation of Ah t and find for Ah t the result 

M ' = -EE ^e-^m,F*D* tli D t2j f(ml m| , m|) 
i=ij=\ 611 

2 2 2 m . 

-EEE g 2 H: 3 {u* kl D* bll - wu^Wi^M.^.^) 
i=i j=i k=i i07r 

2 2 2 m +m + 

i=l j=l k=l 

2 2 4 mo 
"EEE 2^.{a tfc Ai, - l tk D t2j }{ft k D* tli + a t kD* t2i } — ^ / (m 2 o , raf . , m| ) 
i=ij=ifc=i iD7r 

4 4 2 m y om y o 
"EEE 2A?.{ayAi* - 7«A*}{/% Au. + a ti D* t2k }^-^f(ml,m 2 x o, mjo) (39) 

i=l j=l fc=l 

where 

g2m t X 4k 



2mvi/ sin /3 
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(3 tk = eQ t X'* k + X 2 * fc (T 3t - Q t sin 2 9 W ) 

COS U\\r 

,n y> 9 -9l sin2 6w V (Art 
1th = eQ t X lk X 2k (40) 

COS u\\r 
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and where Q t — § and T 3t = \. Similarly for 5h t we find the result 

a f = -EE ^e-^m~ g E tJ D* tu D t2 ,f(ml m|, mf \) 
1=1 j=i 671 

2 2 2 m , 

" E E E 9 Gij{UkiDl u - fea^SsX W^O^/Kj , m? , m|) 
«=i j=i k=i 

2 2 2 777, +777, + 

- E E E ^q,{^^i fe - W% 4, mj+) 

2=1 J = l fc = l 

2 2 4 mo 
- E E E 2£y WAij - itkD t 2j}{l3* tk D* tli + atfcD^}— ^f(m 2 x o, mf ' m] ) 
i=ij=i fe= i 



4 4 2 m Y om y o 

" E E E 2r?.{oyAi* - 7 t ,A2*}{A*A*i fe + "w^l-^/K, mjo, mjo) ( 41 ) 
i=i j=i fc=i iD7r 



4 Correction to the r Lepton Mass 

In this section we discuss the loop corrections to the r lepton mass. The analysis of the 
loop corrections to the r lepton mass are very similar to the analysis of the b quark mass 
with one important distinction. Unlike the correction to the b quark mass where there 
are contributions arising from the gluino, chargino and neutralino exchanges, for the case 
of correction to the r lepton one has contributions arising only from the chargino and 
neutralino exchanges. Keeping this distinction in mind we follow closely the procedure of 
the analysis of b quark mass. Thus we begin by defining the r mass at the Z scale by 

m T {M z ) = h T (M z )-^ cos/3(l + A T ) (42) 
V2 

where A T is loop correction to m T . A T is given by 

,ReAh T n Re8h T , , 
A T = [— tan (3 + —^} (43) 

The remaining loop analysis is very similar to that of the b quark mass and we omit the 

details and give the results below. The analysis of Ah T gives 

2 2 m +m + 

A ^ = " E E 9 2 CijVak T U* 2 * *' f(ml , , m^ + ) 
i=i j=i 

2 Tn v + 
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2 2 4 jjj q 

- E E E 2G* Tij {a Tk D Tlj - lTk D T2j }{(3* Tk D* Tli + a Tk D* 2i }-^f(m 2 x „, m? , m|) 

i=l j=l fc=l iD7r 

4 4 2 m v om y o 

E E E 2r^{a rjJ D rlfc - 7ri A^Hfolk, + q^D^}-^/^, rojo, m 2 o) (44) 



where 



m 



y/2mw cos /3 
gm T X 3k 

2m w cos /3 



/3 Tfc = e g T Xi* fc + -Jl—X'* k (Tz T - Q T sin 2 6 W ) 

COS fvi/ 

Irk = eQ T X lk X 2k (45) 

COS v\\r 



and where 



G Tij _ gM z 1 2 , 



y/2 2cos6»^ 2 



{(-- + sin 2 9 w )D* li D Tlj - sin 2 9\yD* 2i D T 2j } sin /3 



and £V is 

— 7= = — snip (47) 

Similarly for Sh T one gets 

2 2 m x +m x + 

$h T = - E E 9 2 K ij V* l k T U* 2 f(ml T ,m 2 x +,m 2 x +) 



fc=i 



2 2 4 m „ 

- E E E 2//^{a TfcJ D Tlj - 7r*^-H)9;*i>; u + «A}^/Ko,™K) 
i=i j=i k=i io7r 

4 4 2 m y om Y o 

2Ay{a Ti A-i* - 7Ti £> rtJfc }{# i £>; ifc + a Ti D* T2k }^Jlf(ml, m 2 , m 2 o) (48) 



=ii=i fc =i 167f2 



where 



= + si ° 2 <>w)D' TU D Tll - sin 2 e w D' Tli D T2j } cos/3 
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and 



V2 



4 cos 6 



gM z 



•w 



cos (3 



(50) 



5 Discussion of results 



We begin with a discussion of the effects of CP phases on on the SUSY QCD and the 
SUSY electroweak correction to the b quark mass. In Fig. 3 a plot of the b quark mass 
correction as a function of 9^ is given for values of tan (3 ranging from 3 to 50. One finds 
that the correction is very sensitively dependent on 9^ as the value of 9^ affects both 
the sign and the magnitude of the correction. Thus the correction can vary from zero 
to as much as 30% in some regions of the parameter space and can also change its sign 
dependent on the value of 9^. In Fig. 4 a similar plot is given as a function of 9^ for 
different values of £ 3 for the case tan/3 = 50. One finds that Arrib is also a very sensitive 
function of the phase £ 3 and as in the case of 9^ this phase can also change both the 
magnitude as well as the sign of the b quark correction. In Fig. 5 we give a plot similar to 
that of Fig. 4 for the case of tan (3 = 5. The conclusions here are very similar to those for 
Fig. 4 except that the over all magnitude of the correction is typically smaller due to the 
smaller value of tan f3 in this case. In Fig. 6 a plot of the b quark correction as a function 
of (Xa for various values of \A \ is given. One finds again a very strong dependence of the 
b quark correction on both as well as on \A \. This strong dependence is less obvious 
than the dependence on 9^ and £3 and arises from their effects on the stop and sbottom 
masses and via the dependence of the matrices that diagonalize the stop and sbottom 
mass matrices on otA and on \Aq\. In Fig. 7 we give a plot of the b quark correction 
as a function of tan (3 for the three cases of Table 1 of Ref. [20] which satisfy the edm 
constraints for the electron, the neutron and for the H g atom for value of tan f3 = 50. The 
corresponding cases where the phases are all set to zero are also plotted. The plot shows 
that the effect of phases consistent with the edm constraints can produce large effects on 
the b quark mass. 

An analysis similar to the above for the correction to the top quark mass is given in 
figures Fig. 8, Fig. 9, Fig. 10, Fig. 11, and Fig. 12. Thus in Fig. 8 we give a plot of the top 
quark correction as a function of the 9 ^ for all the same parameters as in Fig. 3. As in 
the case of the b quark correction here also one finds that the correction to the top quark 
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mass is very strongly dependent on 9^. However, unlike the case of the b quark correction 
the correction to the top quark mass is typically less than a percent which means that the 
SUSY loop correction to the top quark mass is likely to be at the level of a GeV or so. 
Such an effect could still be significant at the level of precision theoretical analyses. Very 
similar conclusions follow from the analysis of Fig. 9, Fig. 10, Fig. 11, and Fig. 12 which 
parallel the analysis of Fig. 4, Fig. 5, Fig. 6, and Fig. 7. As in the b quark correction here 
also one finds that the SUSY QCD and SUSY electroweak corrections to the top quark 
mass are strongly dependent on the phases £ 3 , otA and on the parameters m , mi and 

Finally, we discuss the numerical size of the correction to the r lepton mass. In Fig. 13 
we give a plot of the correction to the r lepton mass as a function of 9^ for all the same 
parameters as in Fig. 3 and in Fig. 8. As in the case of the correction to the b quark 
mass and to the top quark mass here also one finds that the correction to the r lepton 
mass is very strongly dependent on 9^. However, the numerical size of the correction is 
significantly smaller than the correction to the b quark mass although generally bigger 
percentage wise than the percentage correction to the top quark mass. Again such effects 
could be relevant in precision physics. In Fig. 14, analogous to the case of Fig. 6 for the 
b quark and Fig. 11 for the case of the top quark, a plot of the r lepton correction as a 
function of cxa for various values of \A \ is given. As in the case of Fig. 6 and Fig. 11 
one finds a very strong dependence of the r lepton mass correction on cxa and on \A Q \. 
In Fig. 15, as in the case of Fig. 7 for the b quark mass and in the case of Fig. 12 for 
the top quark mass, we give a plot of the r lepton correction as a function of tan/5 for 
the three cases of Table 1 of Ref.[20] which satisfy the edm constraints for the electron, 
the neutron and for the H g atom for value of tan (3 = 50. The corresponding cases where 
the phases are all set to zero are also plotted. The plot shows that the effect of phases 
consistent with the EDM constraints can produce large effects on the r lepton mass. 
Fig. 15 shows that the loop correction to the r lepton mass can become as large as 5%. 
Thus overall we find that consistently in all cases considered the supersymmetric QCD 
effects and the supersymmetric electroweak effects are strongly dependent on CP phases. 
They can change both the sign as well as the magnitude of the one loop correction. Such 
effects will have important implications on low energy phenomena where supersymmetry 
QCD and supersymmetric electroweak effects enter. An example of this is the analysis of 
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branching ratio for the decay B® d — > l + l~ decay which is strongly affected by the SUSY 
QCD corrections [20, 33] and is also very sensitive to the CP phases. As pointed out in 
Sec.l, Yukawa coupling unification is sensitive to SUSY QCD and SUSY electroweak loop 
correction and thus will also be sensitive to the CP phases. 

6 Conclusion 

In this paper we have computed the effects of CP phases on the b and t quark masses 
and the r lepton mass. In Sec. 2 we discussed the general technique for the computation 
of QCD and electroweak loop corrections to the fermion masses and specifically analysed 
the contributions to the b quark mass arising from the gluino, chargino and neutralino 
exchange contributions. We also checked that our results agree with previous analyses in 
the limit when phases are all set to zero and tan/5 is taken to be large. In Sec. 3 we ex- 
tended the analysis to compute corrections to the top quark mass arising from the gluino, 
chargino and neutralino exchange contributions. In Sec. 4 we extended the analysis to 
compute the SUSY electroweak corrections to the r lepton mass arising from chargino 
and neutralino exchange contributions. These analyses were carried out allowing for the 
most general allowed set of CP phases and without a constraint on tan (3. Thus the 
analysis presented in this paper is valid not only for large tan/5 but also for moderate 
and small values of tan/5. The analysis presented in Sees. 2, 3 and 4 is the first complete 
one loop result of the corrections to the b and t quark masses and to the r lepton mass 
allowing for all allowed phases and valid for all tan/5. In Sec. 5 we gave an exhaustive 
numerical analysis of the size of the SUSY QCD and SUSY electroweak corrections to 
the b and t quark masses and to the r lepton mass. One finds that the numerical size 
of the corrections to the b quark mass can be as large as 30% while the correction to 
the r lepton mass as large as 5%. However, we find that the correction is typically less 
than a 1% for the case of the top quark mass. Further, in all cases analyzed one finds 
that the correction is sharply dependent on the CP phases. Thus one finds that both 
the magnitude as well as the sign may be affected by the presence of the CP phases. 
The loop corrections are of great importance in unified models specifically those involving 
Yukawa coupling unification, e.g., b — r unification for SU(5) and b — t — r unification for 
a class of SO(10) models. More generally these loop corrections will be of import in the 
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study of quark-lepton textures in the presence of CP phases. The SUSY QCD and SUSY 
electroweak corrections to the b quark mass is also of great importance in the analysis of 
low energy phenomena where the b quark mass enters prominently such as phenomena 
involving B mesons. Thus, for example, the decay B® d ^ is strongly affected by the 
SUSY QCD and SUSY electroweak corrections [20, 33] . 
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Figure 3: Plot of the b quark mass correction Amj/mj as a function of 9^ The curves 
in descending order at the point 6^ = it correspond to tan/3 = 3,5,10,20,30,40,50. 
The other input parameters are: m = mi = 200 GeV, £i = .5, £ 2 = -659, £ 3 = .633, 

2 

otA = 1-0, and \Aq\ = 4. All angles here and in succeeding figure captions are in radians. 
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Figure 4: Plot of the b quark mass correction Amj/m;, as a function of 9^ for various 

values of £3. The other input parameters are: m = mi = 200 GeV, £1 = .5, £2 = .659, 

2 

= 1.0, \Aq\ = 4, and tan/3 = 50. The curves in descending order at the point 9^ = 
correspond to £ 3 = 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. 
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Figure 5: Plot of the b quark mass correction Am^/nib as a function of 6^. for various 

values of £3. The other input parameters are: m = mi = 200 GeV, £1 = .5, £2 = -659, 

2 

aA = 1-0, \ A Q \ = 4, and tan/3 = 5. The curves in descending order at the point 6^ = 0.5 
correspond to £ 3 = 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. 
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Figure 6: Plot of the b quark mass correction Am^/m;, as a function of ctA for various 
values of \A \. The other input parameters are: m = mi = 200 GeV, £i = .5, £2 = -659, 
£3 = .633, 6> M = 7r/4, and tan/5 = 50. The curves in descending order at the point 
ctA = 2.5 correspond to \A \ = 1,2,3,4,5,6. 



23 



30 




-50 1 1 1 1 1 1 1 1 ' 

5 15 25 35 45 

tun [3 



Figure 7: Plot of the b quark mass correction Am^/nib as a function of tan/5 for the three 
cases (a), (b) and (c) of Table 1 of Ref.[20] where the EDM constraints for the electron, 
the neutron and H^ 99 atom are satisfied for the case tan (3 = 50. Among the curves in the 
lower half plane corresponding to negative correction to the b quark mass, the long dashed 
curve corresponds to the ) where m = mi = 200 GeV, A = 4, oia — 1, £i = -5, 

£ 2 = -659, £ 3 = .633; the solid curve corresponds to the case (b) where m = mi = 370 
GeV, A = 4, «a = 2, £i = .6, £2 = .653, £3 = .672; the dashed curve corresponds to 
the case (c) where m = mi = 320 GeV, A Q = 3, «a — -85 £1 — -4, £2 = -668, ^3 = .6. 
Ofi = 2.5 in all cases. The three similar curves in the upper half plane are for the three 
cases above when the phases are all set to zero. 



24 




Figure 8: Plot of the top quark mass correction Am t /m t as a function of 9 ^ for values 
of tan/3 ranging from 3 to 50. The other input parameters are: tjiq = mi = 200 GeV, 

2 

£i = .5, £2 = -659, £3 = .633, a^ = 1-0, and \A \ = 4. The curves in ascending order at 
the point 9^ = correspond to tan (3 = 3, 5, 10, 20, 30, 40, 50. 
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Figure 9: Plot of the top quark mass correction Am t /m t in percentage as a function of 0^ 
for various values of £3. The other input parameters are: m = mi = 200 GeV, £1 = .5, 
£2 = .659, ctA = 1-0, \A \ = 4, and tan/3 = 50. The curves in descending order at the 
point 6^ = n correspond to £ 3 = 0, 0.5, 1.0, 1.5, 2.0, 2.5, 3.0. 
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Figure 10: Plot of the top quark mass correction Am t /m t in percentage as a function 



of On for various values of £3. The other input parameters are: mo 



mi 



200 GeV, 



6 = -5, £2 = -659, a Ao = 1.0, \A \ 
the point 6^ = n correspond to £ 3 



4, and tan (3 = 5. The curves in descending order at 
0,0.5,1.0,1.5,2.0,2.5,3.0 



27 



0.75 




O 0.5 1 1 .5 2 2.5 3 

(« u ,(ra<lian) 



Figure 11: Plot of the top quark mass correction Am t /m t in percentage as a function of 
aA for various values of \A \. The other input parameters are: m = mi = 200 GeV, 
£i = .5, ^2 = -659, £3 = .633, 9^ = n and tan/3 = 50. The curves in ascending order at 
the point oca = 2.5 correspond to \A \ = 1, 2, 3, 4, 5, 6 
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Figure 12: Plot of the top quark mass correction Am t /m t in percentage as a function of 
tan/3 for the three cases (a), (b) and (c) of Table 1 of Ref. [20] where the EDM constraints 
for the electron, the neutron and Hjj 39 atom are satisfied for the case tan (3 = 50. The long 
dashed curve with intercept at tan/3 = 50 of 0.36 is for the case (a) where m = mi = 200 
GeV, A = 4, o>a — 1, £i = -5, £2 = -659, £3 = .633. The solid curve with intercept 
at tan/3 = 50 of 0.61 is for the case (b) where mo = mi = 370 GeV, A = 4, o>a = 2, 
£1 = .6, £2 = .653, £3 = .672. The dashed curve with intercept at tan/3 = 50 of 0.49 is for 
the case (c) where m = mi = 320 GeV, A = 3, «a = -8, £1 = .4, £ 2 = .668, £3 = .6. 
# M = 2.5 in all cases. The other three similar curves have all the same parameters as cases 
(a), (b), and (c) above except that phases are set to zero. 
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Figure 13: Plot of the tau lepton mass correction Am T /m T as a function of 9^ for values 
of tan/3 ranging from 3 to 50. The other input parameters are: m = mi = 200 GeV, 

2 

£i = .5, £2 = -659, £3 = .633, a^ = 1-0, and |^4o| = 4. The curves in ascending order at 
the point 9^ = correspond to tan (3 = 50, 40, 30, 20, 10, 5, 3. 
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Figure 14: Plot of the tau lepton mass correction Am t /m t in percentage as a function of 
a.A for various values of \Aq\. The other input parameters are: = mi = 200 GeV, 
£1 = .5, £2 = .659, £3 = .633, 9^ = n and tan/3 = 50. The curves in ascending order at 
the point oca = 2.5 correspond to \A \ = 1, 2, 3, 4, 5, 6 
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Figure 15: Plot of the tau lepton mass correction Am T /m T in percentage as a function of 
tan/3 for the three cases (a), (b) and (c) of Table 1 of Ref. [20] where the EDM constraints 
for the electron, the neutron and Hg atom are satisfied for the case tan/5 = 50. The 
upper dot-dashed curve is for the case (a) where mo = mi = 200 GeV, A = 4, cka = 1, 
£i = .5, £2 = -659, £3 = .633. The upper solid curve is for the case (b) where m = mi = 
370 GeV, A = 4, ctA — 2, £1 = .6, £ 2 = .653, £3 = .672. The upper long dashed curve is 
for the case (c) where m = mi = 320 GeV, A Q = 3, aA = -8, £1 = .4, £ 2 = -668, £3 = .6. 
# M = 2.5 in all cases. The lower three similar curves have all the same parameters as cases 
(a), (b), and (c) above except that phases are set to zero. 
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